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Abstract 

Bounds on the expansion coefficients of the strangeness changing Kir form 
factors were derived recently from analyticity and unitarity, using as input 
suitable correlation functions calculated by perturbative QCD in the Euclidian 
region. We investigate two types of invariant amplitudes and their correspond- 
ing dispersion relations, and show that they lead to similar results for the shape 
of the vector and scalar K13 form factors. 

1 INTRODUCTION 

The semileptonic decays K — > irlvi (known also as decays) represent the gold 
plated channel to extract the element \V US \ of the CKM matrix. The decay rate is 
parametrized in terms of two hadronic form factors that can not be calculated di- 
rectly in perturbative QCD. Various methods, including Chiral Perturbation Theory 
(ChPT), lattice simulations and dispersion relations have been considered recently, 
providing valuable information on the shape of the form factors in the physical region 
and increasing the precision of | V us | extraction (for recent reviews see [U [2] ) ■ 

A fruitful approach to the study of the form factors, proposed some time ago [U 
Hj , is based on the remark that an upper bound on a weighted integral of the modulus 
squared of the form factors along the unitarity cut is sometimes available from 
independent sources. This condition can be exploited through complex analysis, 
leading to constraints on the values at interior points or on the coefficients of the 
Taylor expansion at zero momentum transfer. Mathematically, the problem belongs 
to the standard analytic interpolation theory for the H 2 Hardy class of functions, and 
is referred to as the Meiman problem [5]. For the Kg^ form factors the method was 
applied in [6]-[12]. Recent studies [TO]-[T2] exploit a modified formalism, proposed 
in [13] , which allows one to include also the phase and modulus of the form factors 
along the elastic part of the unitarity cut. 

The starting point of the approach is a dispersion relation that relates an 
invariant amplitude calculated by perturbative QCD in the Euclidian region to 
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the hadronic form factors on the unitarity cut. Since the dispersion relation re- 
quires subtractions, it is convenient to consider derivatives of the invariant functions 
parametrizing the correlator of two strangeness-changing weak currents. The choice 
adopted in [8]-[12] is suitable especially for the form factors of heavy mesons [HI [15]. 
Alternative definitions of the invariant functions and the corresponding dispersion 
relations have been investigated in [6] and a comparison of two definitions in the case 
of the scalar form factor was performed in [7j. The purpose of the present paper is 
to compare the predictions of two choices of correlation functions for both the vec- 
tor and scalar form factors, using the modified version of the formalism, which uses 
information on the phase and modulus on the elastic part of the unitarity cut. Also, 
we use the more precise input available at present, which includes the calculation 
of the correlators to four loops in perturbative QCD P2]-[l9] and the recent precise 
ChPT and lattice calculations [U [2] . 

In section [2] we introduce the form factors and in section [3] we define the 
invariant amplitudes of the QCD correlator of the strangeness changing current and 
the dispersion relations satisfied by them. In section [J] we briefly review the method 
of unitarity bounds and in section [5] we compare the predictions of the two types of 
dispersion relations for the shape of the vector and scalar form factors. 



2 K i3 FORM FACTORS 

Form factors are defined as usual by the matrix element 



{Ap,)\v^\k{ Pk )) = (p^+ PK )v /+(*) + ipK -jv)a. f (t) , (i) 

where = S'j^u is the hadronic strangeness-changing weak current and t = (px — 
Ptt) 2 is the squared momentum transfer. The vector form factor f+(t) represents the 
P-wave projection of the crossed channel matrix element (0| V^|-ff7f), whereas the 
5-wave projection is described by the scalar form factor 

Several low-energy theorems based on flavour and chiral symmetry are available 
for the K(% form factors. At t = 0, where by construction /o(0) = /+(0), 577(3) 
symmetry implies /+(0) = 1. Deviations from this limit are expected to be small 
[20j and have been calculated in chiral perturbation theory and more recently on 
the lattice [HE]. In the case of the scalar form factor, current algebra [21] relates 
the value at the Callan-Treiman point A^tt = — M% to the ratio Fk/F^ of the 
decay constants: 

f (A Kn ) = ^ + A CT . (3) 

To one-loop in ChPT in the isospin limit Act = —3.1 x 10 -3 [22], and the higher 
order corrections are expected to be small [2]. 
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Unitarity and causality imply that the form factors f+(t) and fo(t) are real 
analytic functions in the i-plane cut along the real axis from the elastic threshold 
t+ = (Mk+M w ) 2 to infinity. According to Watson's theorem [23], below the inelastic 
threshold t- m the phase of each form factor is equal (modulo tt) to the phase shift of 
the corresponding partial wave of ttK elastic scattering. Thus, we can write 

f + (t + ie) = |/ + (t)| e *i /a W t+ < t < t in , (4) 

1 /2 

where 5/ (t) is the phase shift of the P-wave of ttK elastic scattering with I = 1/2. 

1/2 

For the scalar form factor a similar relation involving the phase shift 5 (t) of the 

1/2 

5-wave holds. The phases 5 10 (t) are available from the Roy-Steiner analysis of 
Kit scattering [23] and the dispersive studies of the form factors |25|-|27|. Recently, 
information on the modulus of the form factors on the same part of the unitarity 
cut was obtained from the decays r — > Knv T measured by Belle collaboration [28] . 

We mention finally the Taylor expansion used in most experimental analyses 
in the physical region < t < (M K - M w ) 2 of K -> n semileptonic decays [29]-[37]. 
These expansions are written usually for the ratios /+ ; o(i) = /+,o(i)//+(0) as 

f (t) = l + ^t+^±-t 2 + ..., / (t) = l + ^t+^-rt 2 + ... (5) 

Stringent bounds on the dimensionless slopes and curvatures A' and A" have been 
obtained recently with the method of unitarity bounds |10]-[12], In the present paper 
we shall discuss the dependence of these constraints on the dispersion relations for 
the correlation functions used in the method. 



3 QCD CORRELATION FUNCTIONS 

We consider the correlator of two weak currents V^: 

i J d 4 xe^ x (0\T {^(x)K(0)t} |0> = (- M 2 + q^)Ui{q 2 ) + q^U (q 2 ). (6) 

According to perturbative QCD, the invariant amplitudes Hi(q 2 ) and Ho(q 2 ) satisfy 
subtracted dispersion relations. It is convenient to write these relations for the 
functions Xi(Q 2 ) an d Xo(Q 2 ) defined as: 

We used here the variable Q 2 = —q 2 , which is positive on the Euclidian axis. 
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Alternatively, one can write down convergent dispersion relations for the slightly 
different invariant amplitudes Xi{Q 2 ) an d Xo(Q 2 )'- 



dQ 2 ' vr J (t + Q 2 ) 2 ' 

- (n 2, a 2 [Q 4 n (-Q 2 )] 2 [°° t 2 ImU (t) 
Xo{Q ) = , an9W = - / at 



(8Q 2 ) 2 "Wo {t + Q 2 f 

The functions Xi,o(Q ) an d Xi,o(Q 2 ) can De calculated by Operator Product 
Expansion (OPE) and perturbative QCD at large Q 2 > 0. At present, calculations 
to order a 4 are available (see |17j-|19j and references therein). The perturbation 
expansions of Xi(Q 2 ) an d Xo(Q 2 ) are 

Xi(Q 2 ) = * [1 + 0-318 a s -0.062a 2 - 0.162 a? + 0.176 at ...1, (9) 

,2 



Xo(<9 2 ) = ^^g^t 1 + 1-80 a, + 4.65 a 2 + 15.0 a 3 s + 57.4 a 1 



while the expansions of Xi{Q 2 ) and Xo(Q 2 ) rea< l 

Xx(Q 2 ) = _ [1 + 0.318 a s +0.166 a 2 +0.205 ajj +0.504 a 4 s ...], (10) 
Xo(<9 2 ) = 3 ^^2gT M ^ I 1 + 1-167 a, + 1.437 a 2 + 2.495 a 3 s + 5.254 a A s ...]. 

Here a s is the running QCD coupling and m UiS are the running masses of the quarks 
at the scale Q 2 in MS scheme. The higher mass corrections and the condensate 
contributions are negligible. 

We note that the vector amplitude is renormalization scale invariant, therefore 
the perturbative expansion of Xi(Q 2 ) can be obtained from that of Xi(Q 2 ) by taking 
the derivative with respect to Q 2 and using the renormalization group equation for 
the running coupling. On the other hand, in the scalar case the derivatives with 
respect to Q 2 must be calculated from the fixed scale perturbative expansions of IIo, 
setting afterwards the scale /j, 2 = Q 2 , since the scalar amplitude is not renormaliza- 
tion scale invariant. 

The connection with the form factors is provided by unitarity: assuming 
isospin symmetry, the contribution of the Kir states to the positive spectral functions 
Imlli(t) and Imilo(i) leads to the inequalities 

imniw ^ ^ [(t " f+)( r' )]3/2 |/+M|2g(t - (n) 

3M- [(t-t+)(t-M] 1/2 |, ( . M 2n(, f x 



x We are grateful to Dr. N. Chetyrkin for useful comments on this subject. 
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where t± = (M K ±M 7T ) 2 . 

By combining the dispersion relations © and (JSj) , with the unitarity relations 
(jlip . we obtain for each form factor two inequalities of the type 

POD 

/ dt p(t)\F(t)\ 2 < I. (12) 
Jt + 

The function F(t) denotes one of the form factors f+(t) or fo(t), p(t) > is a 
positive semi-definite weight function (depending on the momentum Q 2 ) and / is 
one of correlation functions ([TJ)- (JHJ) , calculated at large Q 2 > by perturbative 
QCD, according to ([9]) and (fTUjh The generic inequality (fT2j) is the starting point 
for deducing analyticity and unitarity constraints on the form factors. 



4 CONSTRAINTS ON THE K £3 FORM FACTORS 

An overview of the formalism was presented recently in [11] . Here we shall give 
a summary of the version proposed in |13j . which exploits an an optimal way the 
information on the phase and modulus for t + < t < i; n , where ti n is the threshold 
for inelastic channels |10|-|12|. 

We introduce the conformal transformation 



2ff\ V "in ~ y/t'm i (13) 

V^in + V^in — t 

which maps the complex i-plane cut for t > t{ n onto the unit disc \z\ < 1 in the 
z = z(t) plane, and define the function g{z) as 
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[z) = F(t(z))[0(z)}- 1 u J (z) W (z). (14) 



Here t(z) is the inverse of the function z = z(t), O(z) is an Omnes function defined 
by 

O(z) = 0(t(z)), 0{t) = exp (I dt^^j , (15) 

in terms of a phase 5{t), which is known from Watson theorem for t < t[ n , and is an 
arbitrary function, sufficiently smooth (i.e. Lipschitz continuous) for t > t; n , and 



„ W = exp | ^EM f d , HOgTI -) . (16, 
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Finally, w(z) appearing in (|14p is an outer function, i.e. a function analytic and 
without zeros in \z\ < 1, whose modulus on the boundary (z = exp(i(9)) is related 
to the weight p(t) appearing in (fT2l) and the Jacobian of the transformation (fT3l) by 



\w(exp(i0))\ 
2^ 



2 



d9 



(17) 
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In the present case the outer functions can be written in a closed analytic form: 
if we use the dispersion relations ((TJ), the outer functions w + (z) and wo(z) for the 
vector and scalar form factors, respectively, have the expressions [6l [8] 



w+(z) 



W (z) 



(1 + i(-Q 2 )) 3 (l - z ~z(t + )f' 2 {l - z z(t_)) 3 / 2 

X (1 - Z 5("Q 2 )) 3 (1 + 5(t+)) 3 /2(l + g(t_))3/2 ■■ 

^v / T^(l + ^) 3/2 



(18) 



V?>{M 2 K - M? 



16^/2^, 

(1 + z(-Q 2 )) 2 (l - z z(*+)) 1/2 (1 - zi(t_))V 2 
(1 - z z(-Q 2 )) 2 (l + 5(t+))V2 (1 + 5(t_))l/2 ' 



Here z is the current variable and is the function defined in f)13|) . 

The alternative outer functions u> +j o(-£) corresponding to the dispersion rela- 
tions ([8]) differ from (|18p by simple factors: 



ti; + (z) 
u) (z) 



o l + z(-Q 2 ) 



z z{ 



-Q 2 ) 



w + (z), 



(19) 



V2 



l-zz(-Q 2 ) 
l+~z(-Q 2 ) 



Wq{z). 



As proven in [13], the function g(z) defined in (|14p is analytic in the unit disc 
\z\ < 1 and the inequality (112p writes as 



where 



1 

2^ 



2tt 



d%(exp(^))| 2 < I', 



[ tin dtp(t)\F(t)\ 2 . 
Jt+ 



(20) 



(21) 



According to the so-called Meiman problem, the relation (|20p leads to constraints 
on the values of g and its derivatives at various points inside the disc \z\ < 1. In the 
general case, consider the value of g{z) and its first K — 1 derivatives at z = 0, and 
the values at A other interior points: 



1 d k g(z) 



k\ dz k 



2=0 



g k , 0<k<K-l; 

e„, z n / 0, 1 < n < A, 



(22) 
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where gk and £ n are given numbers. Then the following inequality holds: 
I 
6 



where 



6v 



6 

z l 


6 

(ziz 2 ) K 


Cat 


l-z 2 x 
(ziz 2 ) k 


1 - ZlZ 2 
{Z2? K 


1 - ZlZ N 

(z 2 z N ) K 


1 - ZlZ 2 


l-zl 


1 - z 2 z N 


(ziz N ) K 


(Z2Z N ) K 


N 


1 - 


1 - Z 2 Z N 


1 - z 2 




1 


A'-l 



>o, 



(23) 



9k zt 



I' 



k=0 



k=0 



(24) 



All the principal minors of the above matrix should also be nonnegative [HE]. An 
equivalent formulation of the condition (|23p was presented in 

The entries of the determinant (I23p are related, by (I14p . to the derivatives 
F&(0), j < K - 1 of F{t) at i = 0, and the values F(t(z n )). Thus, we obtain from 
(12311 constraints on the values of the form factors and their derivatives at t = 0. 



5 RESULTS 

We compare the two types of correlation functions ([7]) and ([8]) by presenting 
the corresponding constraints on the coefficients of the Taylor expansions ([5]), which 
control the shape of the form factors in the semileptonic region. In the numerical 
calculations we took Q 2 = (2GeV) 2 and ti n = (lGeV) 2 , which is a conservative 
choice as discussed in Refs. [ID]- [12]. The perturbative expansions ([9]) and (fT0|) were 
evaluated with a s (2GeV) = 0.308 ± 0.014 [38J. The phases tf[ =1/2 (i) and ffj =1/2 (t) 
for t < ti n are taken from [24 1, while above ti n we used smooth phases approaching 
7r at infinity. As proven in [II] and checked numerically with high precision, the 
results are independent on the parametrization of the phase for t > i; n . Finally, we 
calculated the low energy integrals in (I2ip with the parametrizations of and 
|/o(£)| obtained from the measured rate of the r — > Kttu decay [28]. More details 
on the input and its uncertainties are given in |10]-[12], 

Using the correlation function \\ defined in (J7|), we obtain from (|23p the fol- 
lowing quadratic inequality for the slope X' + and the curvature A'_J_: 

/+(0)[(A'|) 2 - 0.107A' + A'| + 2.18 x Kr 4 A'| + 2.98 x 10- 3 (A' + ) 2 

- 1.49 x 10~ 5 A' + + 4.20 x 10~ 8 ] - 4.67 x 10~ 7 < 0. (25) 
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Using alternatively the correlation function xi denned in (|8|), we obtain 

f% (0)[(A'|) 2 - O.IOA'+A'j + 0.9 x 10~ 4 A'j + 2.59 x 10~ 3 (A' + ) 2 

- 7.44 x 10~ 6 A' + + 2.64 x 1(T 8 ] - 2.99 x 1(T 7 < 0. (26) 

These constraints are shown in Fig. Q] as the interior of two ellipses in the slope- 
curvature plane. For convenience we used the input /+(0) = 0.962 quoted in pQ. 
One may notice that the domains are very similar, especially in their central regions. 




x 10 2 



Figure 1: Allowed domains for the slope and curvature of the vector form factor, us- 
ing as input the correlation function X\{Q 2 ) from ([7]) - solid line, and the correlation 
function Xi(Q 2 ) from (|SJ) - dashed line. 



For the scalar form factor we shall include an additional input using Callan- 
Treiman theorem ([3]). Using the correlation function xo from ([7]), we obtain from 
(f23j) the condition 

/+(0) [(Aq) 2 + 0.25A' Ag + 21.6 x lO^Ag + 15.3 x HT 3 (A' ) 2 

+ 2.68 x 10" 3 Aq + 1.17 x 10~ 4 ] - 10" 3 /+(0)/ct(2.67Ao 

+ 21.53Aq + 0.23) + 1.16 x 10~ 4 /ct - 3-23 x 10" 10 < 0. (27) 

where we denoted fcT = fo(^Kn)- The alternative dispersion relation (JSJ) for the 
correlation function xo leads to 

/+(0) [(Aq') 2 + 0.24A' A' / + 21.1 x KT 3 Ag + 14.4 x KT 3 (A' ) 2 

+ 2.52 x 10~ 3 Aq + 1.12 x 10~ 4 ] - 10- 3 / + (0)/ct(2.53A' 

+ 21.02Aq + 0.22) + 1.11 x 10~ 4 fc T - 2.64 x 10~ 10 < 0. (28) 
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Figure 2: Allowed domains for the slope and curvature of the scalar form factor, using 
as input the correlation function Xo{Q 2 ) from ([7]) - solid line, and the correlation 
function Xo(Q 2 ) from (|HJ) - dashed line. 



The numerical coefficients in these constraints are very close, showing that the two 
types of correlation functions give very similar results. This is illustrated also in 
Fig. [21 where we represent the domains (J2"7|) and (f28j) in the slope-curvature plane 
for the standard input /+(0) = 0.962 and /o(Axtt) = 1.193. 

6 CONCLUSIONS 

The method of unitarity bounds, especially in the improved version proposed 
in [13J . was shown to be a very useful tool in the phenomenological study of the 
form factors [10]-[T2] . In the present paper we investigated the dependence of the 
bounds on the QCD correlation functions and the dispersion relations used as input 
in the formalism. The results show that two different types of correlation functions 
lead to similar constraints for the coefficients that control the shape of the vector 
and scalar form factors at low energies, providing a nice consistency check of the 
formalism. 

Acknowledgements 

I.C. thanks B. Ananthanarayan, S. Ramanan, Gauhar Abbas and I. Sentitemsu 
Imsong for a fruitful collaboration on the subject of form factors. This work 
was supported by CNCSIS in the Program Idei, Contract No. 464/2009, and by 
ANCS in the frame of the project PN 09370102. 



9 



References 

[1] M. Antonelli et al. [FlaviaNet Working Group on Kaon Decays], 
larXiv:1005.2323l [hep-ph]. 

[2] A. Kastner and H. Neufeld, Eur. Phys. J. C 57, 541 (2008). 

[3] S. Okubo, Phys. Rev. D3, 2807 (1971). 

[4] V. Singh and A. K. Raina, Fortsch. Phys. 27, 561 (1979). 

[5] N.N. Meiman, JETP 17, 830 (1963). 

[6] C. Bourrely, B. Machet and E. de Rafael, Nucl. Phys. B189, 157 (1981). 

[7] C. Bourrely and I. Caprini, Nucl. Phys. B722, 149 (2005). 

[8] R.J. Hill, Phys. Rev. D74, 096006 (2006). 

[9] G. Abbas, B. Ananthanarayan, Eur. Phys. J. A41, 7 (2009). 

[10] G. Abbas, B. Ananthanarayan, I. Caprini, I. Sentitemsu Imsong and S. Ra- 
manan, Eur. Phys. J. A44, 175 (2010). 

[11] Gauhar Abbas, B. Ananthanarayan, I. Caprini, I.S. Imsong and S. Ramanan, 
Eur. Phys. J. A45, 389 (2010). 

[12] Gauhar Abbas, B. Ananthanarayan, I. Caprini and I.S. Imsong, Phys. Rev. 
D82, 094018 (2010). 

[13] I. Caprini, Eur. Phys. J. C13, 471 (2000). 

[14] L. Lellouch, Nucl. Phys. B479, 353 (1996). 

[15] C. Bourrely, I. Caprini and L. Lellouch, Phys. Rev. D79, 013008 (2009). 

[16] A. L. Kataev and V.V. Starshenko, Mod. Phys. Lett. A10, 235 (1995). 

[17] PA. Baikov, K.G. Chetyrkin, J.H. Kuhn, Phys. Rev. D67, 074026 (2003). 

[18] P.A. Baikov, K.G. Chetyrkin, J.H. Kuhn, Phys. Rev. Lett. 96, 012003 (2006). 

[19] P.A. Baikov, K.G. Chetyrkin, J.H. Kuhn, Phys. Rev. Lett. 101, 012002 (2008). 

[20] M. Ademollo and R. Gatto, Phys. Rev. Lett. 13, 264 (1964). 

[21] C.G. Callan and S.B. Treiman, Phys. Rev. Lett. 16, 153 (1966). 

[22] J. Gasser and H. Leutwyler, Nucl. Phys. B250, 517 (1985). 



10 



[23] K.M. Watson, Phys. Rev. 95, 228 (1954). 

[24] P. Buettiker, S. Descotes-Genon, B. Moussallam, Eur. Phys. J. C 33, 409 
(2004). 

[25] B. Moussallam, Eur. Phys. J. C 53, 401 (2008). 

[26] B. El-Bennich, A. Furman, R. Kaminski, L. Lesniak, B. Loiseau and B. Mous- 
sallam, Phys. Rev. D79, 094005 (2009). 

[27] V. Bernard, M. Oertel, E. Passemar and J. Stern, Phys. Rev. D80, 034034 
(2009). 

[28] D. Epifanov et al. [Belle Collaboration], Phys. Lett. B654, 65 (2007). 
[29] A. Lai et al. [NA48 Collaboration], Phys. Lett. B604, 1 (2004). 
[30] A. Lai et al. [NA48 Collaboration]. Phys. Lett. B647, 341 (2007). 
[31] O. P. Yushchenko et al, Phys. Lett. B581, 31 (2004). 

[32] T. Alexopoulos et al. [KTeV Collaboration], Phys. Rev. D70, 092007 (2004). 
[33] F. Ambrosino et al. [KLOE Collaboration], Phys. Lett. B636, 166 (2006). 
[34] E. Abouzaid et al. [KTeV collaboration], Phys. Rev. D74, 097101 (2007). 
[35] A. Lai et al. [NA48 Collaboration], Phys. Lett. B647, 341 (2007). 
[36] F. Ambrosino et al. [KLOE Collaboration], JHEP 0712, 105 (2007). 
[37] E. Abouzaid et al. [KTeV collaboration], Phys. Rev. D81, 052001 (2010). 
[38] S. Bethke, Eur. Phys. J. C64, 689 (2009). 



11 



